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Abstract 

Theoretical model for the radiation linewidth in a multi-fluxon state of a long 

Josephson junction is presented. Starting from the perturbed sine-Gordon 

model with the temperature dependent noise term, we develop a collective 

coordinate approach which allows to calculate the finite radiation linewidth 

due to excitation of the internal degrees of freedom in the moving fluxon chain. 

At low fluxon density, the radiation linewidth is expected to be substantially 

larger than that of a lumped Josephson oscillator. With increasing the fluxon 

density, a crossover to a much smaller linewidth corresponding to the lumped 

oscillator limit is predicted. 
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A long Josephson junction is an example of a distributed nonlinear oscillator which 
time-dependent response is associated with its intrinsic spatial dynamics. An externally 
applied magnetic field H penetrates into the junction in form of Josephson fluxons (solitons), 
each of them carrying one magnetic flux quantum $o- As schematically shown in Fig. [I], 
fluxons move across the junction under the influence of the bias current generating an 
electromagnetic radiation at the junction boundary. The frequency / of the radiation is 
given by the Josephson relation / = V/§o, where V is a dc voltage induced by the fluxon 
motion. 

For a lumped (short) Josephson junction, the linewidth 5f of the emitted radiation is 
determined by thermal fluctuations of current passing through the junction. Assuming a 
Nyquist noise spectrum, for the current-biased lumped Josephson tunnel junction the full 
linewidth at half power is given by expression |l],|2| 

5f = Af T .f = -^ r — , (1) 

where k-g is Boltzmann's constant, T is the temperature. The linewidth depends on the 
differential resistance Rd = dV/dl at the junction bias point and the static resistance 
Rs = V/I, where I is the bias current flowing through the junction. 

There are two main regimes which characterize the fluxon motion in long junctions. 
First, a shuttle-like resonant fluxon motion gives rise to zero-field steps (ZFSs) in the dc 
current- volt age (I—V) characteristics of the junction. In this regime fluxons and antifluxons 
undergo reflections from the junction boundaries and the radiation frequency is determined 
by the junction length L and the fluxon velocity v as Jzfs — v/2L . Second, in the high 
magnetic field, the so-called flux-flow regime occurs and is manifested by a flux-flow step 
(FFS) on the I — V curve. In this regime fluxons are created at one boundary of the junction 
and annihilate at the other boundary. The radiation frequency / = v/d& is determined by 
the spacing between moving fluxons d&. In general, for both ZFS and FFS regimes, the 
radiation linewidth 5f should be related to thermal fluctuations of the fluxon velocity v. 
Joergensen et al. |§ obtained a striking general result that, in spite of a different nature of 
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the phase slippage in small and long junctions, the linewidth of the resonant single-fTuxon 
radiation in a long junction is given by the same Eq. (|1|), except for a missing factor of 4 
due to the modified Josephson relation / = V/2$o a t the ZFS. Experiment || showed a 
reasonable agreement with that theory, though some excess linewidth broadening has been 
seen. For the multi-fluxon state which forms FFS it can be argued that internal degrees 
of freedom in the moving fluxon chain may yield a significant contribution in 5f. Here, in 
contrast to the resonant single-fluxon case, local variations of the fluxon spacing d& change 
the radiation frequency. Recent FFS radiation linewidth measurements by Koshelets et al. 
|D showed the scaling of as predicted by Eq. ([!]) but with an effective temperature 
T e Q being by factor of 8 larger than the physical temperature of their experiment. Since 
there was no theory for the radiation linewidth Sf in the flux-flow regime, the reason for the 
excess noise was not resolved. 

In this paper we present a theoretical model for the radiation linewidth in the flux- 
flow regime for long Josephson junctions. We start from the perturbed sine-Gordon model 
with the spatially and temporarily dependent noise current. Using the collective coordinate 
approach we calculate the finite radiation linewidth due to the internal degrees of freedom 
in the moving fluxon chain. The obtained analytical result is evaluated for the relevant 
parameter range of experimentally studied flux-flow oscillators. 

First, we introduce main characteristic quantities for the system. Mean frequency of 
Josephson oscillations for a long junction in the flux-flow regime is 

(fl = ^^W, (2) 

where v is the velocity of the fluxon chain, (H) is the average magnetic field, A = 2Xl + t is 
the effective magnetic thickness of the junction (Al is the London penetration depth, and t 
is the insulator thickness). The quantities essential for the linewidth problem are the mean 
square deviation of the frequency / 
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and reduced r.m.s. linewidth in the flux-flow regime 



^ - 1 ■ (3) 
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In addition to the thermal fluctuations of the fluxon velocity || , there exists an additional 
mechanism of the line broadening in the FFS regime. According to Eq. (|3]), this mechanism 
is directly related to the irregularities of magnetic field distribution in the junction, i.e. 
to the fluxon density fluctuations in the moving fluxon chain. The most general physical 
origin for these fluctuations is the chain's deformation under the influence of thermal noise. 
Therefore, the thermal noise causes two different contributions to the total linewidth in 
the FFS regime. In the following, the dependence of the fluxon chain deformations on the 
velocity and magnetic field and their contribution to the linewidth are calculated. 

The magnetic field distribution in the junction is proportional to the gradient of the phase 
difference (f)(x,t) between superconducting electrodes: H(x,t) = ($o/27rAAj)0a;(x, t), where 
Aj is the Josephson penetration depth. In the following we note the normalized magnetic 
field as h(x,t) = <f) x (x,t). The phase evolution is governed by the perturbed sine-Gordon 
equation (written in the laboratory coordinates) || 

4>xx - 4>tt - sin = 7 + a4> t + n(x, t) (4) 

where the length x and time t are normalized to Aj and to the inverse plasma frequency uj^ 1 , 
respectively. The first term in the right hand side of Eq. (£|) 7 represents the external bias 
current density J normalized to the critical current density J c , a is the damping coefficient, 
and the term n(x, t) represents the thermal fluctuations with the white-noise spatio- 
temporal correlator 

(n(x 1 ,t 1 )n(x 2 , t 2 )) = t=^—S(x 1 - x 2 )6(t 1 - t 2 ) = ar S(x 1 - x 2 )S(t 1 - t 2 ). (5) 

Here the brackets (...) stand for the statistical averaging (over the ensemble), and Eq = 
8hJ c WXj/2e is the rest energy of a fluxon, W is the junction's width and e is the electron 
charge. 



The phase distribution in the flux-flow regime is given by the cnoidal-wave solution 
written in the laboratory coordinates as 

4> = 4>cn(x — vt — £), where (f) cn (z) = tx — 2am(z / h\J 1 — v 2 ). (6) 

Here axn(z) is the elliptic amplitude, v is the velocity of the fluxon chain, and the slowly 
varying shift £ describes the chain's deformation. In the coordinate frame moving at the 



velocity v, x' = (x — vt) / yl — v 2 and t' — (t — vx) / yl — v 2 , the deformation is governed by 
the linear equation J?| 

iff - ix'x' + "6' = p~ lr n(x', f), (7) 

where p = AE(k)/k 2 K(k) is the density of the fluxon chain, the elliptic modulus k is related 
to the chain's period V = 2kK(k) in the moving frame, K(k) and E(k) being the complete 
elliptic integrals of the first and the second kind, respectively. The average velocity of the 
fluxon chain v under the external bias current 7 is given by the Marcus and Imry perturbation 
approach || as 

-1/2 

(8) 
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This relation determines the form of the current-voltage (J — V) characteristics in the flux- 
flow regime. In the laboratory coordinate frame k in Eq.(^|) is determined by given fluxon 



spacing I as a root of the transcendental equation l/yl — v 2 = 2kK(k) ||. 
The local magnetic field in the junction h(x,t) is given by the expression 

h(x, t) = - Vt - £)(1 - £c) ^ [<f>'cn{ x - vt ) - 4>L{ X ~ vt )(] (1 - 6») > ( 9 ) 

where, using the smallness of £, we have Taylor expanded the first multiple in Eq. (0), 
keeping the first two terms. Next we define the mean and mean squared magnetic fields as 

(h) = 0l n (x - vt), (h 2 ) = ^ - fa - fa x ) 2 "j . (10) 

Substituting ( |TUD into Eq. (Q), we obtain the following relation between the reduced r.m.s. 
linewidth and the deformation of the fluxon chain £: 
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Thus, in order to calculate AfpF we need to find the solution of Eq. (|7|) which can be done 
by means of the Fourier transform. First, we write Eq. (|7D in the laboratory reference frame: 

iu - i xx + R - i x + R - i t = p- l n{x, t). (12) 
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Using the Fourier transform t) = f_™ £(g, ui) exp(iqx — uot) dqduj it is straightforward 
to find the solution to Eq. ([12]): 



^ UJ ) = — 2 —( \fi aw/2' ( 13 ) 

q A — uo z — iOi\0J — vq)[ \ — v z ) L > z 

where n(q,uj) is the Fourier transform of the thermal noise n(x,t) = n(q, u) exp(iqx — 
iujt) dqdu which is subject to correlations 

(n(q, u)n(q', u')) = ^-6(q + q')6(u + J). (14) 

Here, we consider only the case of sufficiently dense fluxon chain in the form 

cn ~ (h) (x - vt) - (hy 2 (1 - v 2 )- 1 M(H) (x - vt)), (15) 

which assumes 

(h) 2 (1 - v 2 ) > 1. (16) 

In this case, the coefficient which enters Eq. (|TID is (0c' n /0cn) 2 — W~ 2 (1 — v 2 )~ 2 . Note, that 
according to Eq. (|12"D we have ^/^ 2 ~ a 2 v 2 {l — v 2 )" 1 . It is easy to see that, if we add to 
the chain stiffness condition ([16]) the additional condition 

a 2 (h) 2 v 2 (l-v 2 ) < 1 (17) 



the first term under the square root in Eq. fllT] ) may be neglected as compared to the second 
term. The condition (|I7]) holds in most practically important case of the underdamped 
junction. Since v 2 (l — v 2 ) is always below its maximum value of 1/4, the conditions ( |T6"P 
and fli~7|) can be combined as 
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For the region (ffBI) we thus obtain 



1 « (h) 2 « ± . (18) 



a/ ff = w- 1 (i - « 2 )-y(e 2 (x,t)> , (19) 

where we have inserted the above approximation for (4>c n / 4>' C ri) 2 ■ For the stiff fluxon chain 
its elliptic modulus k and mass density p are given by particularly simple expressions k = 



2/ (h) VT^ 2 and p = (hy (1 - v 2 ). 

Now, what remains is to calculate the r.m.s. value J (£ 2 ). Using Eqs. C[3~3f) and fllil) this 
quantity can be calculated explicitly. However, we encounter a divergence when performing 
integration over small wave numbers q. This divergence is regularized by the fact that, in 
a Josephson junction of the finite length L, the smallest wave number is tt/L. Taking this 
into account we finally obtain 



A fpp = *,ML*f. (20) 

This expression provides the linewidth of a Josephson flux-flow oscillator. It is related to the 
intrinsic mechanism of the fluxon chain deformation under the influence of thermal noise. 

The calculation of another noise contribution A fa due to the thermal motion of the 
fluxon chain as a whole (similar to that of Ref. 0) give the result identical to that of 
Eq. (|l]). Thus, the total reduced linewidth of the Josephson flux-flow oscillator is given by 
the expression 

Afe = Af T + Af FF . (21) 

Let us discuss the physical meaning of the result (|2T|). First, the term fc^T / 'Eq is typically 
small, of the order of magnitude 1CT 4 — 10 ~ 5 . Therefore, a comparison of Eq. ([1]) with 
Eq. ( PUD shows that, due to the square-root dependence of AfpF on ksT/Eo, the linewidth 
in the flux-flow regime can be essentially broader as compared to the resonant flux motion 
regime. Moreover, the flux-flow linewidth scales with the junction length as vL, which 
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means additional broadening for longer junctions. Physically, this dependence is related to 
a formal divergence of the chain's deformation in the infinitely long junction. 

According to Eq. (p0|), the linewidth is rapidly reduced with the increase of the magnetic 
field h which makes the fluxon chain more stiff. In sufficiently high magnetic fields h > h* 
the contribution determined by Eq. ( ^) may become smaller than that due to thermal 
fluctuations of the chain velocity given by Eq. (0). In this case the crossover to the standard 
mechanism discussed in || takes place and the total linewidth is determined by Eq. (fj). 

The crossover magnetic field h* depends essentially on the dimensionless ratios ksT/aEo, 
L/Xj and on the fluxon chain velocity v. We note that Eqs.(|I]) and (|20|) predict quite different 
dependencies on v. Namely, A fx decreases with v and becomes very small in the relativistic 
regime v — > 1. In contrast to that, Afpp increases with v. 

In order to compare Afr and AJff it is useful to estimate them for typical parameters of 
practical Josephson tunnel junctions. Figure | shows the calculated reduced linewidths given 
by Eqs. (|IJ) and (|20| ) as functions of the fluxon velocity for two typical sets of experimental 
parameters H^j. In order to determine Ro and R$ in Eq. ([I]), the form of the current- 
voltage (/ — V) characteristics for the fluxon spacing I = 2ir/h was calculated using Marcus 
and Imry perturbation formula (|5p which with given parameters was independently checked 
by numerical simulations for periodic boundary conditions. In both cases we took the dissi- 
pation coefficient a = 0.02 and the normalized magnetic field h = 4 (which is approximately 
twice the first critical field of the fluxon penetration into the junction) at T = 4.2K. In case 
of low-J^ junctions || (Jc = 200A/cm 2 ) we used L = 200/im, W = 10/im, Aj = 35/im, 
and in the other case of high- Jc junctions [||] (Jc = 8000A/cm 2 ) the values L = 450/im, 
W = 3/im, A j = 4/iin were used. In Fig. 2 one can see that in both cases the effect of the 
spatially-dependent fluxon chain noise contribution Afpp given by Eq. ( |20"D is dominating, 
in particular at high fluxon velocities. At low velocities the reduced linewidth Af? (|1|) is 
formally diverging due to finite value of 5 fx at / = 0, while the intrinsic flux-flow linewidth 
AfpF saturates at the finite level according to Eq. (^). In both cases the rest fluxon energy 
ksT/Eo is of the order of 10 _5 J. For the used magnetic field h = 4 the stiff fluxon chain 
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assumption ([16]) strictly holds only for v < v s — 0.6. Since at v > v s the fluxon chain can 
only become more soft, Eq. (^) gives the lowest limit for the expected Af FF . 

We note, however, that in experiment the flux-flow oscillators are often operated close to 
a resonant regime. In such cases the junction frequency spectrum is modulated by character- 
istic eigenfrequencies which determine the oscillation amplitudes (P~5|) . In the resonant states 
the linewidth can be expected to be narrowed by cavity resonances (known as Fiske steps) 
in the junction. Thus, the non-resonant flux-flow frequency linewidth given by Eq. (^0|) can 
be taken as the upper limit for the real flux-flow oscillator of the finite length. 

Finally, we like to mention that additional mechanism of the flux-flow oscillator linewidth 
broadening may arise from technological inhomogeneities in the junction. Such a contribu- 
tion is nonuniversal and depends strongly on the particular type of disorder (local imper- 
fections in the tunnel barrier, the precision of the photolithography-defined junction width, 
etc.), and therefore this mechanism was not discussed in the present paper. This issue will 
be the subject of further investigation. 
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FIGURES 

FIG. 1. Schematic cross-section of a long Josephson junction in the flux- flow state. The elec- 
tromagnetic radiation emitted at the junction boundary depends on the variations of the fluxon 
spacing due to temporaly and spatially dependent noise current through the junction. 



FIG. 2. Reduced flux-flow oscillations linewidth 5f = Af/f given by Eq. (20) (solid lines) and 



Eq. (|H) (dashed lines) as a function of the fluxon chain velocity v. Thick and thin lines correspond 
to two typical parameter sets accounting for Refs. || and §, respectively. 
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